GatekeepR: an R Shiny application for the identification of nodes with high dynamic impact in Boolean networks.

MOTIVATION: Boolean networks can serve as straightforward models for understanding processes such as gene regulation, and employing logical rules. These rules can either be derived from existing literature or by data-driven approaches. However, in the context of large networks, the exhaustive search for intervention targets becomes challenging due to the exponential expansion of a Boolean network's state space and the multitude of potential target candidates, along with their various combinations. Instead, we can employ the logical rules and resultant interaction graph as a means to identify targets of specific interest within larger-scale models. This approach not only facilitates the screening process but also serves as a preliminary filtering step, enabling the focused investigation of candidates that hold promise for more profound dynamic analysis. However, applying this method requires a working knowledge of R, thus restricting the range of potential users. We, therefore, aim to provide an application that makes this method accessible to a broader scientific community. RESULTS: Here, we introduce GatekeepR, a graphical, web-based R Shiny application that enables scientists to screen Boolean network models for possible intervention targets whose perturbation is likely to have a large impact on the system's dynamics. This application does not require a local installation or knowledge of R and provides the suggested targets along with additional network information and visualizations in an intuitive, easy-to-use manner. The Supplementary Material describes the underlying method for identifying these nodes along with an example application in a network modeling pancreatic cancer. AVAILABILITY AND IMPLEMENTATION: https://www.github.com/sysbio-bioinf/GatekeepR https://abel.informatik.uni-ulm.de/shiny/GatekeepR/.

A protocol for the use of cloud-based quantum computers for logical network analysis of biological systems.

Boolean networks are commonly used in systems biology to dynamically model gene regulatory interactions. Here, we present a protocol for implementing Boolean network dynamics as quantum circuits. We describe steps for accessing cloud-based quantum processing units offered by IBM and IonQ and downloading and parsing logic for gene regulatory networks. We then detail procedures for performing simulations of quantum circuits on local devices and visualizing measurement results. For complete details on the use and execution of this protocol, please refer to Weidner et al.1.

A systems biology approach to define mechanisms, phenotypes, and drivers in PanNETs with a personalized perspective.

Pancreatic neuroendocrine tumors (PanNETs) are a rare tumor entity with largely unpredictable progression and increasing incidence in developed countries. Molecular pathways involved in PanNETs development are still not elucidated, and specific biomarkers are missing. Moreover, the heterogeneity of PanNETs makes their treatment challenging and most approved targeted therapeutic options for PanNETs lack objective responses. Here, we applied a systems biology approach integrating dynamic modeling strategies, foreign classifier tailored approaches, and patient expression profiles to predict PanNETs progression as well as resistance mechanisms to clinically approved treatments such as the mammalian target of rapamycin complex 1 (mTORC1) inhibitors. We set up a model able to represent frequently reported PanNETs drivers in patient cohorts, such as Menin-1 (MEN1), Death domain associated protein (DAXX), Tuberous Sclerosis (TSC), as well as wild-type tumors. Model-based simulations suggested drivers of cancer progression as both first and second hits after MEN1 loss. In addition, we could predict the benefit of mTORC1 inhibitors on differentially mutated cohorts and hypothesize resistance mechanisms. Our approach sheds light on a more personalized prediction and treatment of PanNET mutant phenotypes.

Leveraging quantum computing for dynamic analyses of logical networks in systems biology.

The dynamics of cellular mechanisms can be investigated through the analysis of networks. One of the simplest but most popular modeling strategies involves logic-based models. However, these models still face exponential growth in simulation complexity compared with a linear increase in nodes. We transfer this modeling approach to quantum computing and use the upcoming technique in the field to simulate the resulting networks. Leveraging logic modeling in quantum computing has many benefits, including complexity reduction and quantum algorithms for systems biology tasks. To showcase the applicability of our approach to systems biology tasks, we implemented a model of mammalian cortical development. Here, we applied a quantum algorithm to estimate the tendency of the model to reach particular stable conditions and further revert dynamics. Results from two actual quantum processing units and a noisy simulator are presented, and current technical challenges are discussed.

CANTATA-prediction of missing links in Boolean networks using genetic programming.

MOTIVATION: Biological processes are complex systems with distinct behaviour. Despite the growing amount of available data, knowledge is sparse and often insufficient to investigate the complex regulatory behaviour of these systems. Moreover, different cellular phenotypes are possible under varying conditions. Mathematical models attempt to unravel these mechanisms by investigating the dynamics of regulatory networks. Therefore, a major challenge is to combine regulations and phenotypical information as well as the underlying mechanisms. To predict regulatory links in these models, we established an approach called CANTATA to support the integration of information into regulatory networks and retrieve potential underlying regulations. This is achieved by optimizing both static and dynamic properties of these networks. RESULTS: Initial results show that the algorithm predicts missing interactions by recapitulating the known phenotypes while preserving the original topology and optimizing the robustness of the model. The resulting models allow for hypothesizing about the biological impact of certain regulatory dependencies. AVAILABILITY AND IMPLEMENTATION: Source code of the application, example files and results are available at https://github.com/sysbio-bioinf/Cantata. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Identification of dynamic driver sets controlling phenotypical landscapes.

Controlling phenotypical landscapes is of vital interest to modern biology. This task becomes highly demanding because cellular decisions involve complex networks engaging in crosstalk interactions. Previous work on control theory indicates that small sets of compounds can control single phenotypes. However, a dynamic approach is missing to determine the drivers of the whole network dynamics. By analyzing 35 biologically motivated Boolean networks, we developed a method to identify small sets of compounds sufficient to decide on the entire phenotypical landscape. These compounds do not strictly prefer highly related compounds and show a smaller impact on the stability of the attractor landscape. The dynamic driver sets include many intervention targets and cellular reprogramming drivers in human networks. Finally, by using a new comprehensive model of colorectal cancer, we provide a complete workflow on how to implement our approach to shift from in silico to in vitro guided experiments.

Reconstructing Boolean network ensembles from single-cell data for unraveling dynamics in the aging of human hematopoietic stem cells.

Regulatory dependencies in molecular networks are the basis of dynamic behaviors affecting the phenotypical landscape. With the advance of high throughput technologies, the detail of omics data has arrived at the single-cell level. Nevertheless, new strategies are required to reconstruct regulatory networks based on populations of single-cell data. Here, we present a new approach to generate populations of gene regulatory networks from single-cell RNA-sequencing (scRNA-seq) data. Our approach exploits the heterogeneity of single-cell populations to generate pseudo-timepoints. This allows for the first time to uncouple network reconstruction from a direct dependency on time series measurements. The generated time series are then fed to a combined reconstruction algorithm. The latter allows a fast and efficient reconstruction of ensembles of gene regulatory networks. Since this approach does not require knowledge on time-related trajectories, it allows us to model heterogeneous processes such as aging. Applying the approach to the aging-associated NF-κB signaling pathway-based scRNA-seq data of human hematopoietic stem cells (HSCs), we were able to reconstruct eight ensembles, and evaluate their dynamic behavior. Moreover, we propose a strategy to evaluate the resulting attractor patterns. Interaction graph-based features and dynamic investigations of our model ensembles provide a new perspective on the heterogeneity and mechanisms related to human HSCs aging.

Capturing dynamic relevance in Boolean networks using graph theoretical measures.

MOTIVATION: Interaction graphs are able to describe regulatory dependencies between compounds without capturing dynamics. In contrast, mathematical models that are based on interaction graphs allow to investigate the dynamics of biological systems. However, since dynamic complexity of these models grows exponentially with their size, exhaustive analyses of the dynamics and consequently screening all possible interventions eventually becomes infeasible. Thus, we designed an approach to identify dynamically relevant compounds based on the static network topology. RESULTS: Here, we present a method only based on static properties to identify dynamically influencing nodes. Coupling vertex betweenness and determinative power, we could capture relevant nodes for changing dynamics with an accuracy of 75% in a set of 35 published logical models. Further analyses of the selected compounds' connectivity unravelled a new class of not highly connected nodes with high impact on the networks' dynamics, which we call gatekeepers. We validated our method's working concept on logical models, which can be readily scaled up to complex interaction networks, where dynamic analyses are not even feasible. AVAILABILITY AND IMPLEMENTATION: Code is freely available at https://github.com/sysbio-bioinf/BNStatic. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.